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1 Introduction 

Consider the initial boundary value problem for the KdV equation on the 
first quadrant 

Ut = Uxxx ~ QuUx, X > 0, t > 0, (1) 

u\x^o = a, Uxxlx^o = b, (2) 

u\t=o = uo{x), wo(a:)|^^+oo ^ 0. (3) 

The initial value uo{x) is supposed to be consistent with the boundary con- 
dition (0) at the corner point: uq{0) = a, woa;a;(0) = 6, the boundary data a 
and b are chosen to be constants. Besides uo{x) should be a smooth func- 
tion decreasing rapidly enough. Remind that the problem (|I])-(@) is correctly 
posed and uniquely soluble (see, for instance, []T|). 

The problem (|I])-(0) admits infinitely many integrals of motion. The first 
three of them are of the form 



/ udx = (3a — b)t + const, 

f°°{ul + 2u^)dx = -{3a^ - b)H + const, 

/ i^xx ~ ^^ ^xx + 5w + 4(6 — 3a )u )dx = const. 

When the parameter 3a^ — 6 is different from zero then evidently the first two 
integrals of motion contain the explicit ^-dependence. It can be proved (see 
[]2j]) that the problem ([I|)-(0) is consistent with an infinite number of higher 
symmetries and, because of this reason, the KdV equation (|I]) admits a large 
class of explicit algebra-geometric solutions satisfying the boundary condition 



see [^). All these facts allow one to hope that the problem (|I])-(0) can be 
integrated by means of the inverse scattering method. The scattering matrix 
s{^,t) of the associated hnear operator y" = {u{x,t) — X)y, x > 0, defined 
in the standard way evaluates in time t by means of the following system of 
equations 

St = 4^4 [s, CTsJ + [Uj^ai — (72 H — la^Js (4) 

containing unknown Wa;(0, t) and given w(0, t) = a and Uxx{0, t) = b (here and 
below aj are the Pauli matrices). The equation (||) shows that the nonlocal 
change of variables from u{x^t) to s{^,t) doesn't lead to any separation of 
variables in the problem (|T|)-(0). Note that such kind obstacles always arise 
when the initial boundary value problem is studied for integrable equations 
(see, for instance, [^, 0] where an alternative approach to the problem is 
discussed). 

The system (g) is underdetermined, both the coefficient and the solution 
are not given (that is why the system is really nonlinear). But it is worth 
to note that in addition to the system there is an extra condition: for all 
values of t the solution s(^,t) has to preserve its analytical properties, i.e. 
has to belong to the class of scattering matrices. The latter allows one in 
some particular cases to simplify essentially the system (^ and ever reduce 
it to a problem of factorisation. Let us do first a linear change of dependent 
variables by setting s(^, t) = T{^)S{^, t) exp{4:i^^(Tst). Then the new variable 



S{^,t) solves the equation 

St = {iiyas + w^cri)^, (5) 



where 



and 



no = I ^. .^1,9 = 4^' + 2a, 



v{0 = ^16^6 + 4^2(5 _ 3^2) ^ 2ab - 4a\ (6) 

Suppose the function P = P(^) to define an involution of the curve (JOj) 
such that z/(^) = z/(P(^)) for all ^. Then evidently the identity 

S-\P{0,t)S{U) = S-\P{O,0)S{^,0) (7) 

holds, which may be considered as the first integral of the system (0). The 
problem arises how to solve the nonlinear equation ((T]) for t > 0. Below 
in the paper we discuss the particular case when the additional constraint 
b = 3a^, a < is imposed and however we act in a slightly different way 
(the homogeneous case a = b = has been studied in [0). Here we find 
the solution of the t-equation at the straight-line a: = by gluing up the 
eigenfunctions of the two operators from the Lax pair and then solving the 
Riemann problem of linear conjugation for piecewise analytic functions on 
a six sheeted Riemann surface. The main observation we use below is as 
the X- and ^/-equations have not only a common fundamental solution but 
a common solution of the scattering problem for all values of the spectral 
parameter except some contour on the complex plane. 



2 The direct scattering problem 

It is well known that the KdV equation is the consistency condition of the 
following two systems of ordinary linear differential equations: 

Y, = UY, (8) 

Yt = VY, (9) 

where the coefficient matrices are 

/on / tt^ -4A-2w 

\u — X 0/ V^xa; ~ (4A + 2w)(w — A) —Ux 

Due to the boundary conditions (@) the equation (g) takes the form 

Yt={ ]y (10) 

V6-(4A + 2a)(a-A) -u^ J 

along the border a: = 0. In other words the linear systems (g), (0), and (pXI|) 
give the Lax representation of the initial boundary value problem (|I]), (g), 
and (0). 

Define matrix- valued solutions of the auxiliary system (g), satisfying the 
following asymptotic for all real values of the parameter ^ = vA 

Y,{x,t,O^To{Oe'^^^\ x^+oo, (11) 

Y2{x,t,O^To{^), x^O, (12) 

/I 1 \ 
where To(^) = . The columns i/jk{x,t,^) and (j)k{x^t,^), k = 1,2 

of the matrices Yi = [i/ji, 0i) and Y2 = (02, '4'2) are known to have analytical 



continuations with respect to the parameter ^ from the real axis onto the 
complex plane. Moreover, Y2{x,t^^) is an entire analytical function and its 
columns have the following asymptotic in the corresponding half-planes when 
1^1 goes to infinity 



^2 



{x,t,0 ^ e-^ f M (l + 0(r')), Im^ < 0, 



(13) 



Mx, t, ^ e-^^^ i-^^^^^ 0(r')), Im^ > 0. (14) 

In general the columns of the other solution Yi are defined only on the 
half-planes Im^ < and Im^ > 0, respectively, where they have similar 
asymptotic for ^ ^ oo 

0i(x, t, ^ e— ^ [ ^ ] (1 + 0(r')), Im^ < 0, (15) 



^i(x,i,0 ^ e^^^ fM (1 + 0(r')), Im^ > 0. 



(16) 



The scattering matrix of the system (@) given on the half-line is defined as 
a ratio of two fundamental solutions: s(^) = Yq^^Yi. It is easy to show that 
entries sn and S21 of the matrix s are analytic in the domain Im^ > and 
similarly the entries su and S22 ~ in the domain Im^ < 0. 

Let us compound matrices and i/j by taking the vectors (J)i and ipi as 
their columns multiplied by scalar factors as follows 

0(x,O = (02,0lS2l')(x,Oe-^"^"^ (17) 

ij{x,0 = (^isn\^2)(x,0e-^"^"^. (18) 
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By definition ttie functions 0(a:,^) and ip{x^^) are analytic in lower and 
upper half-planes, respectively. At the real axis Im^ = they are related to 
each other by conjugation condition: 

c^{x,0=i^{^,Or{x,0- (19) 

The X dependence of the conjugation matrix r(a:, ^) is given for x > by: 

r{x, = e*^^'^3r(0, ^e'^^^^^^ (20) 

The asymptotic representations above yield 

0(x, -^ ToiO for Im^ < 0, <e ^ oo and (21) 

^{^,O^To{0 for Im^> 0,^^00. (22) 

The eigenfunctions (|r7|), (|TB|) are degenerate at the point ^ = 0, since the 
scattering matrix has generally a pole there. Evidently, 

0(x,O)Q=O, ilj{x,0)r]=0. (23) 

In the special case when s(^) is bounded at ^ = the condition ([23| ) is not 
valid. 

The conjugation matrix in (|T^) is expressed explicitly through entries of 
the scattering matrix: 

/ 1 Si2S^2^e2*^^\ 

\-S2iSue ''""^ SnS22 J 
A non-degenerate matrix-valued solution Y{x, ^) of the system (g) is called 
solution of the scattering problem on the half-line a: > at the point ^, Im^ 7^ 



0, if the function Y{x,^)e~^^^^^ is bounded uniformly for all values of a: > 
0. Clearly functions Y{x,^) = 0(a:,^)e*^^'^3 ^^d Y{x,^) = ip{x,^)e^^^'^^ are 
solutions of the scattering problem at the corresponding half-planes except 

the points ^ where S22{0 — 0? ^^ I™^ < ^^^ ^uiO ~ 0? i^ I™^ > 0- 

Remark. The solution of the scattering problem on a half-line is not 

unique. It is defined up to the right side matrix-valued multiplier upper 

triangular when Im^ > and lower triangular when Im^ < 0. 

We will use this freedom when finding solutions of the similar scattering 

problem for the second operator which is responsible for the t-evolution. 
Let us give explicit representations of the matrices (j) and il) at the end 

x = Q 

^0-^(00(0,0= "f hro"'(OV^(0'0=L^ J- (25) 

3 Method of gluing up the eigenfunctions 

In the previous section we constructed eigenfunctions (i.e. solutions of the 
scattering problem) for the x-equation (g). According to our conjecture there 
should be a common eigenfunction for both x- and t-equations (§|) and (pXI|) at 
the corner point (x = 0, t = 0) defined for all values of the spectral parameter. 
But how to find such a solution? The matter is that the scattering problem 
for these two equations are naturally formulated on two different spectral 
planes. Really, after the linear change of variables Y{t,^) = T{^)Y{t,^) (pXI|) 



takes a form of the well-known Zakharov-Shabat system (cf. (|5|)) 

Yt={iiyas + u,{0,t)ai)Y, (26) 

where the function 



jy = ^16^6 + 2a3 (27) 

is a particular case of (0) under the constraint b = 3a^. 

To be able to compare functions on ^ with those on u it is necessary 
to consider two Riemann surfaces F and H appearing in a natural way such 
that the function (|27j) establishes a one-to-one correspondence between them. 
The first one is a two-sheeted cover of the ^-plane glued by its six cuts such 
that two cuts go along the real axis four others are done to make surface 
unvariant under rotation by angle | (see Pic.l). Denote through Fi and F2 
the sheets of the surface. Require that ^ G Fi if the asymptotic representation 
iy{^) = 4^^ + 0(1) is vahd for ^ ^ 00 and < arg^ < |. Remind that the 
continuous spectrum of the equation (g) coincides with the real axis Im^ = 0. 
By rotating the axis around the origin by angles | and ^, i.e. applying the 
involution of the curve, one gets two more directed lines. The conjugation 
contour L, consisting of these three lines taken on both sheets, divides the 
surface into twelve sectors. Enumerate them in the following way. Denote 
through /i 1 and /2,i those lying in the first quadrant on the sheets Fi and F2 
respectively. Denote /2,i and /2,2 the next two sectors that lie partly in the 
first and partly in the second quadrants on the corresponding sheets and so 
on. 



The function ^ = ^{i/) inverse to ([27|) is evidently single-valued on the 
six-sheeted cover of the complex plane H with sheets Hj, j = 1, ...,6 glued 
by two cuts lying on the imaginary axis and going to infinity from the branch 



points iy± = ±iy2|aP (see Pic. 2). The sheets are enumerated in such a way 
that the images Ikj = i^{hj) are as follows 

/^ ^ = {i/ e Hk, Rez/ > 0} U {z/ G Hk+i, Rez/ < 0} n {Imu > 0}, 
/2^^ = {ly e Hk, Reiy < 0} U {z/ e Hk+i, Rez/ > 0} n {Imiy < 0}. 

The sheets are such that when crossing the upper cut from the left to the 
right one passes from Hj to Hj+i, where Hj = Hi and Hq = Hq. 

Denote through Ykj{t, v) = ^kj{t, v)e~^^^ ^^ solutions of the equation ([2B|) 
where the multipliers ^fcj(t, z/) are regular for large enough u from Ikj. Ac- 
tually, these multipliers just are to be found. Normalize them by fixing their 
triangular matrix structure 

(^ \ / ^ ^ \ 

, *fc,2m+l(0,Z/)= . (28) 

Comparison of these eigenfunctions with those of the x-equation gives rise to 
the so-called gluing up equations: 

^(0, 04,(0 = T{0^k,M ^)^kjH (29) 

forj = 1,2,3; /c= 1,2; 

^(0, 04,(0 = r(O*fc,(0, iy)Akj{iy) (30) 

for J =4,5,6; k=l,2. 
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The left side of these equations gives an expression of arbitrary solution of the 
scattering problem for the x-equation while the right side gives that of the 
other one. Here the right side multipliers Skj{^) are upper (lower) triangular 
matrix-valued functions defined on Ikj for j = 1,2,3 {j = 4, 5, 6, respectively) 
and /S.kj{v) are upper (lower) triangular for j = 2,4,6 (j = 1,3.5) matrices 
defined on I^j. Thus two multipliers 5kj and A^j are of one and the same 
matrix structure for only four choices of subindices {k,j): (1,2), (2, 2), (5,1), 
(5,2). In these cases the equations ([2^ ) and ( |3D|) are effectively solved by 
reducing to the problem of triangular factorisation of matrices. For example, 
when /c = 1, j = 2 one has 

The rest of functions are defined by imposing the following involution 

^kMO) = ^k,j+2{i^H)) = ^kM^i^'O). (31) 

which is completely consistent with the equations ([2^ ), (^) as well as the 
involutions 

*fcj-(-i^(-0) = ^l*mn(^(0)^l (32) 

and 



^^kji-H-O) = ^^mniHO), (33) 

where the line over a letter means the usual complex conjugation. 

Having all of the functions "^kjit, t^) at the point t = one can find the 
conjugation matrices of the Riemann-Hilbert-Carleman problem on the Rie- 
mann surface H with the contour of jumps L = v{L) associated with the 
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^-equation from the Lax pair: 

^fcj(t, ly) = "^mnit, iy)Rkjmn{t, l^) ■ (34) 

The conjugation matrices exphcitly depend on t: 

Rkjmnit, V) = e-^'^^'^jkn^niO, v)e^'^'^' (35) 

By solving the Riemann problem one finds the eigenfunctions ^kj{t-, ^) for all 
t > and recovers the linear system (|2B|), i.e. its unknown coefficient Wa;(0, t). 
Afterwards the equation for the time evolution of the scattering matrix (|5D 
becomes linear. On the other hand side the scattering matrix s(^,^) can 
directly be expressed in terms of the functions s(0, ^) and "^kj{t^ ^)- And the 
last step consisting of recovering the potential u{x^ t) for a given s(^, t) called 
the inverse scattering problem is studied very much (see, for instance, [jTj]). 

To conclude we remark that since the spectrum of t-equation for x = 
lies on a two genus Riemann surface then it is natural to expect that Wa;(0, t) 
has finite-gap behaviour at infinity. 
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Pic.l. The double sheeted Riemann surface F. 
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Pic. 2. The six sheeted Riemann surface H 
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